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rG ' Abstract 

ix . The asymmetric tetrahedron is composed by all edges of tetrahedron represented by Ising interaction except 

one, which has a Heisenberg type interaction. This asymmetric tetrahedron is arranged connecting a vertex which 
edges are only Ising type interaction to another vertex with same structure of another tetrahedron. The process 
is replicated and this kind of lattice we call the asymmetric Ising-Heisenberg chain. We have studied the ground 
state phase diagram for this kind of models. Particularly we consider two situations in the Heisenberg-type 
interaction, (i) The asymmetric tetrahedral spin(l/2,l/2) Ising-XYZ chain, and (ii) the asymmetric tetrahedral 
spin-(l/2,l) Ising-XXZ chain, where we have found a rich phase diagram and a number of multicritical points. 
Additionally we have also studied their thermodynamics properties and the correlation function, using the dec- 
orated transformation. We have mapped the asymmetric tetrahedral Ising-Heisenberg chain in an effective Ising 
chain, and we have also concluded that it is possible to evaluate the partition function including a longitudinal 
external magnetic field. 
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1 Introduction 

>, 

^sO , Low dimensional systems based on magnetic material have attracted considerable attention lately in a number 

0^ i of subjects such as condensed matter physics, material science and inorganic chemistry. In these particular areas 

quantum ferrimagnetic chains (QFC) were discussed, due to that they exhibit a relevant combination of ferromag- 
netic (F) and antiferromagnetic (AF) states. Experimental synthetization of the compound CuG-Clpy^fN^fTll 
with Clpy indicating Chloropyridine had been investigated. This compound could be mapped into a spin-1/2 
tetramer chain with F-F-AF-AF bond alternation[2j. Recently diamond type chain structures have been intensively 
investigated theoretically and experimentally (3J. The natural candidates to describes these kinds of materials are 
the quantum anisotropic Heisenberg model or even Ising type models. Certainly the rigorous mapping of those 
compounds into Heisenberg type models could become very complex systems which usually have non exact solu- 
tion. However some particular case of models could become exactly solvable such as the Ising-Heisenberg model 
considered by Jascur and Strecka[4], a more detailed discussion also considered by Canova et al.|5|. The method 
used to solve this kind of models is the historical Fisher's work[6) of decorated transformation method, proposed 
in 50's decade, the improvement of the method is discussed in reference |7|. Many other quasi-unidimensional 
Ising type models were solved using this methodsl8l[9). 

Recently, theoretical investigation of strongly geometrical frustration materials[10] have been performed, par- 
ticularly focused on the diamond chain structure, using several numerical approaches fill [l2l . These theoret- 
ical results could enhance the other experimental realizations provided by the polymeric compounds such as 
Cu20S04|13) and M3(OH)2 (with M=Ni, Co, Mn)|HHT5|. Other quasi-unidimensional Heisenberg models were 
studied using numerical results such as discussed in references fl2l[T6ll and some analytical series expansion also 
has been performed[17J for similar systems. 

This work aim is to present the frustrated properties of the asymmetric tetrahedron Ising Heisenberg (ATIH) 
chain, this models can be solved exactly mapping for an effective Ising chain with spins 1/2 or 1 using the method 
presented by Fisher[6|. This work is organized as follow, in section 2, we present the ATIH chain, considering 
Heisenberg interaction with spin 1/2 and 1. In section 3, we discuss the phase diagram at zero temperature 
shown a rich phase diagram and several critical points, for the spin-(l/2,l/2) Ising-XYZ chain and spin-(l,l/2) 
Ising-XXZ chain. In section 4, we discuss the thermodynamics properties for Ising-Heisenberg chain with Ising 
spin s=l/2 or 1 and Heisenberg interaction spin S'=l/2 or 1, we also considered the correlation function using the 
decorated transformation method[6J. Finally in section 5, we present our conclusions. 
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2 The model 

The asymmetric tetrahedral Ising-Heisenberg (ATIH) chain is composed by all edges of tetrahedron (dashed line 
Fig.l) represented by Ising type interaction except one which is represented by Heisenberg type interaction (solid 
line in Fig. 1), which can be viewed also as Ising-Heisenberg diamond chain. To obtain the ATIH chain we coupled 
the vertex composed only by the Ising type interaction edge, we call from now on just as Ising interaction vertex, 
and connecting to another Ising interaction vertex of other tetrahedron. On the other hand the asymmetric edge of 
the tetrahedron are represented by Heisenberg type interaction, that we simply call as Heisenberg interaction edge 
(interaction between sites a, bin Fig.l). 




Figure 1: The schematic representation of the coupled asymmetric tetrahedral Ising-Heisenberg. All edges 
(dashed line) of tetrahedral are represented by the Ising interaction vertex except one (solid line) which is rep- 
resented by Heisenberg interaction edge. 

The schematic representation of the ATIH model is given in Fig.l. The Hamiltonian for the ATIH chain dis- 
cussed above could be written by the following expression 
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with s ci being the spin of Ising interaction vertex, where J is the interaction parameter. The last term included in 
eq.lQ} correspond to the Heisenberg interaction edge, which is given by 



TjXYZ _ J qX rri , T qV qV , T qZ qZ 
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being S\ the spin matrices with v = x,y,z and J v are their interaction parameter among sites a and b. We can also 
include the longitudinal external magnetic field in the Hamiltonian (l}, which read as 



H„ 



E 



/'o 



\8c,'< 



>c,i+l) 



^ \Sa,i + £>b 



(3) 



here ho is an external magnetic field acting on spin s c , whereas h is an external magnetic field acting on S* and 
S§. Note that we are considering different external magnetic fields because we assume that the gyromagnetic 
factor acting on s c could be different from that gyromagnetic factor acting on S* and S§, which we report as 
ho = fjh with g being their relative gyromagnetic factor. The eq. <[TJ) is symmetric Hamiltonian in relation to the 
exchange s ci ^ s c ,i+i an d S* i ^ S^. On the other hand we note that the Hamiltonian ||lj also has an internal 
spin symmetry H{s c ^ s c , i+l ) = H(-s C}i , s C}i+1 ). 



2.1 The XYZ interaction edge with spin-1/2 

To perform the partial summation over decorated site, we need to diagonalize the XYZ interaction edge Hamilto- 
nian I0, to evaluate this, we introduce the notations J+ = J x + J y and J_ = J x — J y . For the spin S = 1/2, we 
obtain the diagonalized Hamiltonian 

a i+1 = diagfAi 1) ,Af,AL 2) ,A« 



(4) 

with diag() we represents diagonal elements of the Hamiltonian (j}, and conveniently we use the notation for 
simplicity s c = s c ,i and s' c = s c ,i+i- Thus the eigenvalues are given by 
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with 

a = a(s c , s' c ) =J (s c z + s' c z ) + h, (6) 

7 = 7 ( Sc , s' c ) =J Sc z s' c z + ^ (s c z + s' c z ) , (7) 

here a(s c , s' c ) and ~f(s c , s' c ) are dependents of the spins s c and s' c , thus we write from now on just as a and 7 
respectively. Each eigenvalues are given by eq.101 which depend of s c and s' c , then we have 16 eigenvalues. 

After diagonalize the Hamiltonian function we get the corresponding set of eigenvectors also, thus the new 
basis are given by 

\v[ + \s c , s' c )) =_L=( ei | + +> + I - -)), \v{-\s c , s' c )) =—L=(e 2 \ + +> + I - -)), (8) 

V 1 + ej V 1 + e 2 

l*4 +) ) =7^ (I+-) + h+» , h4 _) ) =^= (- l+-> + h+» , (9) 

where the factors e\, e 2 depends of spins s c and s' c which are given by the following values 



/ 16a 2 + j2+4a -Jl6a 2 + J 2 _ + Aa 
ei = ei(s c , s' c ) = ■* , e 2 = e 2 (s c , s' c ) = , (10) 

the normalized eigenvectors \v[ (s c ,s' c )) and |t; 2 ) belonging to the eigenvalues A^ and X 2 given by (0. 
We also remark that the eigenvectors \v. 2 ) are independents of the spins s c and s' c or we can say it is four-fold 
degenerated. 

At this point we would like to comment about the relations of the ei(s c , s' c ) and e 2 {s Cl s' c ) factors. From (TlOb , it 
is not difficult to note that these factors transform one eigenvector of eq.© into the other one, when we exchange 
the following values, 

^\v[ + \ Sc ,s' c ))^\v[-\s c ,s' c )), (11) 

this property could be valid in more general situation, even in the presence of an external magnetic field. However, 
it is important to point out that the eigenvalues of eq.10 does not have a similar transformation. The eigenvalues 
A^ will always be at higher energy level than the eigenfunction A_ . In the case of null magnetic field (h = 
0, ho = 0) we have these coefficients given by 

ei (±, T ) = -e 2 (±,T) = l, (12) 
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ei (±,±) = ^ , e 2 (±,±) = — X . (13) 

These properties will be useful when we discuss the phase diagram properties. 

2.2 The XXZ interaction edge with spin-1 

Now we can consider XXZ interaction among sites a and b (see fig.l), with spin 5=1 case. We diagonalize the 
Hamiltonian analogous to the previous case. After diagonalizing the Hamiltonian depends only of s c and s' c and 
read as 

H iii+1 =diag(A« Af ,Af ,\<?,\M,XW,\?\\™,\V) , (14) 

with diag() we represents diagonal elements of the Hamiltonian |(l]l, whereas the eigenvalues are given by 
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where a and 7 already was defined by the eqs.© and 10. 

Now for a complete analyze of the ATIH model, we shall turn our attention to study the XX Z interaction 
edge for the decorated spin-1. In this situation we have nine eigenvector after diagonalizing the Hamiltonian, 
corresponding to the eigenvalues l|15lll9|l , thus the normalized eigenvectors are 

|n+> =|1,1), \u~)=\-l,-l) 7 (20) 

| U +)=-±=(|l,0) + |0,l)), |«2>=-^(|0,-l> + 1-1.0)), (21) 
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where f± and f 2 are respectively given by 

J z + y /J?+8J* J z - y/J* + 8g 

h = 57 1 h = 57 . ( 2t >) 

at this point, we have a similar situation as was pointed out in eq. (ffl , i.e., these factors transform as 

Jz—>-Jz, Jx—>-Jx then \u%) -»■ luj), (26) 

again this conclusion is valid even when an external magnetic field is included. 

After rewritten the Heisenberg interaction edge Hamiltonian in the diagonal form we are able to discuss the 
phase diagram for whole quasi-unidimensional chain for both cases XYZ interaction edge with spin-1 /2 and XXZ 
interaction edge with spin-1. 

3 The Phase diagrams of the ATIH chain 

3.1 The asymmetric tetrahedral spin-(l/2,l/2) Ising-XYZ chain 

To study the phase diagram of the asymmetric tetrahedral spin-(l /2,1 /2) Ising-XYZ chain we use the diagonalized 
version of the Hamiltonian presented in the previous section 2.1. We would like to note that for the present 
considered model, we have a set of sixteen different state vectors, nevertheless, when the translation and the 
global spin inversion symmetry are taking into account, we obtain that only eight state vectors have different 
phases. 

As we mentioned above some state vectors are physically equivalents, for example, the state vectors \v[ (+,+)) 
and \v[ (— , — )) corresponding to the same state when we consider the global spin inversion. Once the eigenval- 
ues satisfy the relations Aj_ (s c , s' c ) = Aj_ (— s c , —s' c ) and Aj_ (s c , s' c ) = a±'(—s c , ~ s ' c )r 'which restrict the eigen- 
vectors to only eight relevant energy states. So we obtain the following eight states vectors for the asymmetric 



tetrahedral spin-(l/2,l/2) Ising-XYZ chain 
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where the first element inside the product corresponds to the Ising interaction taking two possible values (±1), 
and the next element represent the XYZ interaction edge considered in the previous section. All products are carried 
out over all spin sites. In these relations the single Ising site magnetization mo is given for the spin-1/2 and mi 
is the single Heisenberg magnetization for the a, b sites (Heisenberg interaction edge). The first two states j27[l and 
J28t are new states and arise when we consider the total Ising-XY Z case and their magnetization mi depends of 
the parameters J, J_ and h. These states we will call in general "Quantum ferromagnetic" (QFO) states of type 
I and II, respectively, which are not degenerated and from l(8]l it is possible to see that the probability of the spin 
alignment, defined by the functions ei(s c , s' c ) and e2(s c , s' c ) are not equivalents for the up and down orientation. 
Worth remark that this "Quantum ferromagnetic" state could become Quantum Ferrimagnetic state when mi 
or m-2 is negative. The energy states of (|29) and IpO) also have the structure of the first two states but now the 
probability of the up and down orientation are equivalent, these states are non-degenerated, we call also these 
state as QFO states of type III and IV respectively. The states <(3TT) - ((34b display four frustrated (FRU) states of types 
I, II, III and IV respectively, and are non-degenerated. For these states we use an extended modified unitary cell 
necessary for the identification of equivalent vector states. 

In Fig. 2, we show the ground-state phase diagram for the system in the absence of the magnetic field (h = ho = 
0). In Fig. 2a let us to consider the following re-parametrization for the interaction parameters, without loosing 
any physical generality, 

J = sin (a;), J z = sin(y), J + = 4cos(y), J_ = 4cos(x). (35) 

Here we represent the interaction parameter by means of two new parameters x and y. For the first two states 
(eq. l27l and 1281 we obtain that the Heisenberg interaction edge magnetization takes the value, mi = sin(x)/2, thus 
for x = 0, 7r we have mi = 0, and for x = it/2 we obtain mi = 1/2. 

Certainly, by using this new parameters we restrict the values of interaction parameters as follows: \J\ ^ 1, 
\J Z \ ^ 1 and | J± | ^ 4. In this limited region we have competing interaction parameters, leading to several ground 
state energies. Out of this region there are no new phases. 

It is possible to see that in this phase diagram, we have five tricritical points where converge three states and 
two four-critical points where converge four states. In the last case we have that these points are located in the 
vertical line x = ir, for this position the value obtained from |(35) is J = 0, J z = 0.8, i.e., in this region the system 
has pure Heisenberg interaction for a, b sites, and in general, this region is defined by the line where a continuous 
phase transitions occurs. We are able to calculate the other interaction parameters corresponding to these two 
critical points (x, y), thus we have that for the points (tt, 0.9273) and (ir, 2.2143) the interaction parameters takes 
the values J x = —1.7, J y = 2.3 and J x = —3.2, J y = 0.8, respectively. In principle we must plot in Fig. 2a for 
the parameters x and y in the interval [0, 2ir}. But in order to highlight the rich region in the interval [0, tt], we 
considered it only up to 37r/2 the x parameter, since we have the continuation of the states QFOm, QFOu and 



QFOiv- On the other hand, in the y axis we have also repeated states (FRUi, FRUn, FRUiv and QFOn) which 
are also present in the interval considered in Fig. 2a. We will use the same considerations in the next figures. 

In the region where the axis takes the values x < it, we have that the interaction parameter J > 0, i.e. it is 
positive defined, and for the region x > it we obtain a negative value for the parameter J < 0. In general negative 
values for interaction parameters favors the anti-parallel alignment for the spin ordered system. Thus define the 
type of phase ordered state whatever is FRU or QFO sates (left and right side at the axis x — it in the diagram). 
At this stage it is necessary to point out that the phase state QFOn appears in both sectors (grey sector of Fig.2a), 
basically this occurs by the fact that the constant e2(s c , s' c ) is different from zero. As will be mentioned below, if 
the constant e2(s c , s' c ) comes close to zero, the region QFOn of the left hand side disappears. 

The other five tricritical points can be found by a simple substitution, for example the point (ir/2, tt/2) is where 
the phase states FRU IH , FRU IV and QFO n coexist. Other tricritical points are (0.9273, 1.4679), (0.9273, 1.6735), 
(2.2143, 1.4679), (2.2143, 1.6735). 
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Figure 2: Schematic representation of the phase diagram, (a) For the ATIH chain using the new parameters x and 
y defined by (55) . (b) For the ATIH chain using the new parameters x and y defined by i 



If we impose the condition J x = J y , i.e. for the XX Z model, the eight states reduce to seven, six of them will 
appear in the phase diagram. This is possible because in the limit (J_ —> 0) the constants ei(a c , s' c ) and e 2 (s c , s' c ) 
become 1 and respectively, doing the states FRUj —> FRUn equivalent. 

We give in Fig.2(b), the ground-state energy for regions resulting from the substitution J_ = and where six 

phases for the ground-states energy are shown. Itiseasy to see that the states vectors become \v[ '(s Cl s' c )) — > | + +) 
and \v[ (s Cl s' c )) — > | ). For convenient reasons we use the following realization 



J = -sin(ar), J x = 2cos(y), J z = -sin(y), 



(36) 



and obtain now a ferromagnetic (FM) and ferrimagnetic (FRI) states. It is possible to see that the x = -k axis 
still present in these phase configuration. On this line three critical points are shown where four ground-states 
converge. The phase diagram displayed in Fig.2, shows several states due to the presence of the Js Ci jS Ci j+i inter- 
action, i.e. the j(s c , s' c ) term (with ho = 0) given by (7) allows us to consider this interaction. In reference |5] the 
term Js c ^s c ^ + i was considered null. If we put j(s c , s' c ) = and at the same time impose the equality of some of 
the interaction parameters, we obtain a very simple phase diagram. This is graphically displayed in Fig.3. The 
Fig. 3(a) shows that in the case j(s c , s' c ) — 0, J x — J y , some states are degenerated, so we identify the QFOiy 
and FRU iv as having the same energy. The same occurs with the QFOm and FRUm states. The FM and FRI 
states also appear in this diagram. The Fig. 3b shows that for the case j(s c> s' c ) = 0, J x = J y = J z we obtain the 
more simple configuration with three phases, the FM, FRI and FRU states. Only one three-critical point where 
these states converges is present in (ir,ir/2). At this stage we would like to remark that a very similar ground 
state configuration was obtained in |5), however, the realizations used in our work (35) and (36) are different from 
those used in ||5j, because we did not include the external magnetic field. 
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Figure 3: Schematic representation of the phase diagram, for ATIH chain, where we consider 7 = and J_ = 0. 
In (a) we display a ferromagnetic FM and ferrimagnetic FRI states, however the quantum ferromagnetic QFO 
and the frustrated FRU phases are also present. In (b), we fixed the condition J x = J y = J z obtaining only three 
states phase energy. The FM, FRI and FRU sectors converges into one tricritical point (it, ir/2). 



3.2 The asymmetric tetrahedral spin-(l/2,l) Ising-XXZ chain 



To obtain all states we need to consider the coupling of the decorated vector states l|20H24b with the Ising interaction 
vertex with spin- 1/2, this enable us to write down the total vector states of the system. We will restrict the condi- 
tions over the interaction parameters without losing the generality, using the following values for the interaction 
parameters 

J = —J z = sin(i-), J x = 2sin(y), (37) 

in the following lines we also restrict the system to the case when the external magnetic field is absent (ho = 0, 
h = 0). From all possible twenty-four ground-states energy, which can be obtained from the system only fifteen 
eigenvalues have different values. In this situation the ground-state eigenvectors for the asymmetric spin-(l/2,l) 



Ising-XXZ chain, that would appear in the phase diagrams are given by 
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where the first element inside the product corresponds to the Ising interaction taking two possible values (±1), 
and the next element represents the XXZ interaction edge considered in the previous section 2.2. The eq. ((38t rep- 
resents a ferromagnetic (FM) state, eq. 1 1391 indicate a ferrimagnetic (FRI) state, whereas eqs. l|40ll -ll42ll corresponds 
to Quantum ferromagnetic (QFO) state of type I, II and III respectively. We also have two types of quantum fer- 
rimagnetic (QFI) states given by eq. l[4"3l and ^44), and finally eq.l l45)l represents a frustrated (FRU) state. All these 
states given by J38H45b are displayed in the phase diagram presented in Fig. 4(a). It is remarkable that the line 
x = 7T divides the QFO states from the other ones. Therefore in this case we have a seven-critical point in (tt, it) 
where all states converge. 

In order order to display the first vector state (|38l l, which not present in fig. 4a, we change the restriction over 
the interaction parameters and fixed as follow 



J = J z = sm(cc), J x = 2sin(j/), 



(46) 



the phase diagram is given in the Fig.4(b), where the FM state appears, the other states were already presented 
in the Fig.4(a). 

Finally it is quite interesting to mention that the state given by eq. (|38l l also appear when the next nearest 
interaction parameter considered is null, such as considered in reference [5] i.e. when 7 = 0, and the restriction is 
extended over the decorated interaction parameters, for example, if we put 



J = - sin(a;), J x = J y = J z = 4sin(y), 7 = 0, 





(47) 



we obtain the more simple phase diagram with four states as displayed in the Fig.5. The FM and the FRI states 
are separated by the x = ix line and they end up in the QFOju and FRUi. Similar phase diagram was presented 
in the work J5) with zero external magnetic field and zero next-nearest interaction parameter. In the work |3 
they also consider an external magnetic field but with zero next nearest interaction parameter where a similar 
phase diagrams were obtained, as displayed in the Fig. 4 and 5(a). So we conclude that the inclusion of the (s c , s c ) 
interaction enable us to investigate a rich number of states even in the absence of an external magnetic field. Both 
restrictions l|46H47l l give the value for the ground-state energy Eq = 0. 
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Figure 4: Phase diagram for the ATIH with spin-(l/2,l) chain. In (a) we have obtained it, in more general case, 
seven phase diagrams, using the parameters given by eq. (|37t . In (b) we have obtained it, restricting the interaction 
parameters to J = J z = sin(ai), J x = J y = 2 sin(y), given by eq.(| 



4 The ATIH chain thermodynamics 

Thermodynamics properties could be studied using the known decorated transformation spin proposed in refer- 
ence HHZ|. Let us write the partition function as follow 
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where N is the number of decorated bounds, whereas the T^{s\ stands for the trace of the central decorated 
system or Heisenberg interaction edge (in our case), while by w(s c , s' c ) we represent the Boltzmann weight. One 
should notice that the transformation (|48l l is rather general, since it is valid for arbitrary spin values contained in 
the decorated site. The set {s c } represents the Ising interaction vertex and would take any spin value too. In this 
section we study the Ising-Heisenberg chain where the Ising interaction vertex take the spin values S = 1/2 or 1. 
The ATIH chain partition function can be expressed as 



z = rz , 



(49) 



where Zq is the partition function of the effective Ising chain with arbitrary spin-5, whereas / means a constant 
for the effective Ising chain. 

When decorated Heisenberg interaction edge sites are occupied by spin-1/2 or 1, it is necessary to perform the 
partial trace over all those decorated spin sites. 

4.1 The Ising interaction vertex with spin-1/2 

In order to map the ATIH chain into an effective Ising chain, let us consider the Ising interaction vertex with spin- 
1/2. For this case the associated Boltzmann weight function w(s c , s' c ) has the form 



w(s c , s' c ) 



e A + +c A + 



(50) 



where A± and \±' are given by eq.10 for S = 1/2 (spin of decorated sites) whereas the associated Boltzmann 
weight for spin S — 1 could be obtained using the eqs. JT5T >- (fT9b , then the ATIH chain model considered here is 
written by 



w(s c ,s' c ) 



2e -/ 3 (T+I J *) cosh ((3\^Jl6a 2 + J 2 _) + 2e-^-3 J =) cosh((3\j + ), S = 1/2 

c -/37( e /3J, + 2C 2 J ; cosh (&Jj* + 2J2) + 4 cosh(f J + ) cosh(/3a) + 2c-' 3 ' 7 - cosh(/32a)), S 



(51) 
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Figure 5: Schematic representation, for the connectors with spin-1/2 and decorated Heisenberg spin-1, of the 
phase diagram for the decorated XX Z spin-1 model showing the presence of the FAI state. We fixed the param- 
eter interaction J x = J y = J z = 4sin(y), J = l/5sin(x) and 7 = 0, obtaining only three states. The FAI, FRI are 
non-degenerated while the QFO and FRU are two-fold degenerated. These sectors converge into one tricritical 
point (ir,ir). 

The effective Ising chain partition function is represented by their Boltzmann weight function w(s c , s' c ) which 
read as 

w (s c , s' c ) = / exp {-/3 (Ks c s' c + B(s c + s' c ))} , (52) 

Using decorated transformation, we obtain the new parameters for the effective Ising chain, 






~/3A"=41n 



1 



w{\^) 2 



2 ' 2 



2 \w( — ,—) 



(53) 



where the new effective parameters of Ising chain can be expressed as a function of the parameter of the original 
Hamiltonian. Thus / is just a constant, whereas K means a coupling parameter and finally B corresponds to the 
external magnetic field. 

The expression for the partition function of ATIH chain results in 



N 



Z = f N Z = f N ^ JJ e -0(Ks c s' c +B(s e +s' c )) 



(54) 



Using the eq. ((54t we are able to map the asymmetric tetrahedral spin-(l/2, S) Ising-Heisenberg chain, into an 
effective spin-1/2 Ising chain. Where the spin S of Heisenberg interaction edge could be 1/2 or 1. 



4.2 The Ising interaction vertex with spin-1 

Another case that we consider will be the Ising interaction vertex with spin-1. Thus the asymmetric tetrahedral 
spin(l,S) Ising-Heisenberg chain, will be mapped into an effective spin-1 chain. Similar to the previous case we 
obtain the following Boltzmann weight expressed as follow 



w(s c ,s' c ) = fexp{-p(K lSc s' c + B(s c + s' c )+D(s c 2 + s c 2 ) +E( 
where K\, K-i, B, D and E are the parameter to be determined. 



3 r 2 s' c + S c s' c 2 )+K 2 S c 2 S 



c 2 )}, (55) 



The new parameters of eq.([55]) can be expressed using the associated Boltzmann weight, which is written as 

f=w(0,0), (56) 



1 ln / w(l,l,)w(-l,-l)w(0 , 0) 2 
4 n l w(l,0) 2 w(-l,0f 



-PK 2 =- A In ( ^^'/ „ 2 , ,2 • (61) 



Similar to the previous case / means just a constant in the new effective Hamiltonian, while K\ being the 
coupling parameter, whereas B corresponds to the external magnetic field, the parameter D represents the single- 
ion anisotropy, whereas E correspond to the interaction of quadratic and linear interaction among the nearest spin 
and finally K 2 are the parameter of the biquadratic interaction. 

In this case if we consider a null magnetic field, we have w(l, 1) = w(— 1, — 1) and w(l, 0) = w(— 1, 0). Under 
this condition the eqs. (|58] l and eq. (f60b leads to B = and E = 0, respectively, then the Boltzmann weight function 
w(s c , s' c ) defined by eq.d55l), reduce the following relation 

w(s c , s' c ) = /exp ( - /3(K lSc s' c + D(s c 2 + s' c 2 ) + K 2 s c 2 s' c 2 )) , (62) 

Finally we have concluded that, our mapping of the asymmetric tetrahedral spin-(l,S) Ising-Heisenberg chain 
can be expressed as an effective spin-1 Ising chain. Where as before the spin S of Heisenberg interaction edge could 
be 1/2 or 1. 

4.3 The ATIH chain correlation functions 

We can notice that the partition function of the ATIH chain obtained above, by mapping into Ising spin chain is 
limited. We cannot obtain directly the correlation function because the mapped Ising chain does not depend of 
the decorated spin. Then we can use the method presented by Fisher|6), where the correlation function for the 
ATIH chain can be obtain using the decoration transformation in a similar way as was performed the partition 
function, assuming we known the correlation function of the effective spin-1 Ising chain. 
Using the definition given in reference |6), we have 

(s iSklSk2 ...) = l y, Y. s > s ^--- e ~ PH > < 63 ) 

here 5, represents the decorated spin at site i, and s k are any spins of the systems along the chain, it could be 
either decorated spins S or undecorated spins s c . With Z as the total partition function of the system and H the 
total Hamiltonian. It is possible to split the above relation in two parts, one of this being independent on the spin 
Si and the other one containing the Si dependence. Thus we write down the total Hamiltonian |Q} as 

H = H (Si, 8 c ,i, s c ,2, ■•■) + H n (s c ,i, s c ,2, ■••) = H + H n , (64) 

where the Hq contains the dependence of the spins (Si, s c ,i, s c ,2, ■■■) an d H n contains only combinations of the 
spins (s c .i, s c ,2: ■■■)■ Actually, in the case when [Hq, H n ] = 0, we have 

(S l s kl s k2 ...) = -^ y^ s kl s k2 ...e^ (3Hn n(s c ^s c .i +1 ), (65) 

and as was pointed out by Fisher[6|, it is possible to prove that 0(s Cj jS C! , + i) can be represented as 

fi(*o,<*c,<+i)=$^S i e H9flb . (66) 



We give some values for the correlation functions of the asymmetric tetrahedral Ising-Heisenberg model in the 
case when the Ising interaction vertex with spins are equal to s c = 1/2 or s c = 1 and the decorated spin could be 
equal to 1/2 or 1. In the case of spin s c = 1/2, it is possible to obtain an equivalent form[6] for the right hand side 
of (gill, 

to(s c> i,s Cl j) ={qo + Qo,i(sc,i + s c ,j) +qijsc.is c .j) ^e' 131 * . (67) 

Si 

As an example, let us apply to evaluate the following correlation, with arbitrary sites i and j, instead of 
performing only among next nearest sites(5|, considering S* could be spin 1/2 or 1, thus the correlation read as 

( S a.i s cj) = Qo(s c ,j) + Q0,l(( S c,iSc,j) + (s c ,i+lS c ,j)) + 3l,l(s c ,iS c ,i+lS c ,j) (68) 

where the coefficients q's can be obtained solving the system equation d65t and i68l , from where we verify their 
solution is given as a derivative of the parameters obtained in (f53l> respect to magnetic field h, which read as 

1 d w ldB a ldK «m 

qo = ~2pdh lnf > qifi= 2-dh' and 91 ' 1 = 8^' (69) 

then it is possible to write as a combination of the correlation function of the effective Ising chain with up to 
three-body spin correlations, this correlation is given explicitly as follow 



\^c,i^c,i-\-l'^c,j 



Sc,j) ={sc) 2 + (1 - (*c> a )e (i - J)/4 (70) 



> =( Sc ) 3 + ( Sc )(l-(s c ) 2 )(e( i -^+e- 1 /«(l + e«- i )/«)) (71) 

where £ is the correlation length of the effective(standard) Ising chain (18| . 

The auto correlation function of the ((S*) 2 ) when we consider the Ising interaction vertex s c = 1/2, simply 
becomes a constant equal to 1/4, but when we consider spin-1, this expression becomes non-trivial and we can 
obtain analogous to previous correlation discussed, then the general expression can be write as 

({S^) 2 ) =Qa + 2Qi, (sc,i) + Qi.i (s c ,i, aa+i) (72) 

onces again the coefficients Q's can be obtained using the following relation 

(73) 
(74) 
(75) 

with 

M{s c>i ,s c>i+1 ) = ^ + w ^ , (76) 

whereas w(s c ,i, s c ,i+i) already was defined in eq.([5TT). Other correlation functions like (S^S^), are null when 
v 7^ v' . Whereas for v = v' , the correlation function can be obtained directly using the derivatives of the free 
energy instead of using the previous iterative method such as the one performed by Canova[5|, particularly we 
show the following nearest correlation, 

(S^SQ = -j^ffi, with u = {x,y,z}. (78) 

We remark that above conclusion are also valid even for spin s c = 1. Similar analysis could be performed to 
obtain the correlation function when s c = 1, in this case the fi(s Ct i, s c ,i+i) function will be defined as 

tt(s Ci ,s c ,j) ={qa + qo,i(s c ,i + S c ,j) + qi,iS c ,iS c ,j + qi,2(s c ,iS 2 cd + s 2 ci s cJ ) + q ,2(s 2 cl + s 2 cj ) 

+?2, 2 < i < j }e-^ ) (79) 
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2/3 2 dh 2 
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using this relation, together with the correlation function of the effective Ising chain with spin-1, we can obtain 
other correlations functions for s c = 1, using the same recipes above. 

Alternatively we can obtain this kind of correlation using the direct transfer matrix formalism such as per- 
formed in reference [19J, meanwhile the advantage of this method could be a non-iterative calculation. 



5 conclusions 

The phase diagrams of the asymmetric tetrahedral Ising-Heisenberg (ATIH) chain were studied for the case when 
the Ising interaction vertex is spin-1/2. Firstly we considered the XYZ interaction edge with spin-1/2, and null 
external magnetic field (h = ho = 0). The diagrams displayed in Fig. 2a, have shown seven states appearing in 
the model, with five critical transition points (x, y) having three phase state converging and two critical transition 
points, where four phase states converge. These states have shown their quantum nature for the XYZ interaction 
edge (decorated sites), for example, for the vector states l l2"7t and l |2"5t we see that up and down orientation in the 
decorated sites have different probabilities defined by the factors e\(s c , s' c ) and e-2(s c , s' c ). We also have analyzed 
the particular case when J_ = 0, constructing the phase diagram and showing that FM and FRI states appear 
(Fig. 2b). Other situations have been also studied, particularly, the case where 7 = 0, J x = J y and for the more 
simple configuration 7 = 0, J x = J y = J z/ (Fig. 3). Secondly, when the XXZ interaction edge (decorated) spin is 
1, we have also obtained a rich phase diagram for the ground-state energy, even when an external magnetic field 
and the next nearest interaction are absent. These phase diagrams are shown in the figures 4 and 5. The x = ir line 
appears to be the limit of the QFO and the other states, so for the critical point (tt, it) the energy takes the value 
Eq = 0. Some particular case of our results obtained have been compared with those obtained in reference 0. 

We have noticed that using a decorated Ising model mapping transformation, initially given by Fisher (6), the 
calculation of the partition function for the ATIH chain is reduced to a closed expression of Ising spin chain. We 
have considered some particular cases to discuss the thermodynamic properties, such as Ising interaction vertex 
with spin 1/2 and 1, whereas the interaction edge could be XYZ with spin-1/2 and XXZ with spin-1 respectively. 
The results for correlation function are presented generally for the situation when we have the Ising interaction 
vertex with spin-1/2 or 1 and the Heisenberg interaction edge with spins 1/2 or 1. We have observed that some 
correlation function could be obtained using the derivative related to some parameter instead of using the dec- 
orated transformation used in reference |5), we also considered a long range correlation function where we used 
the decorated transformation to obtain the result. 
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